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Abstract of corresponding document: US2005075274 

A method for determining the interaction of one or more chemical active compounds with organisms. The 
method comprises selecting at least one compound and analyzing its physiko-chemical and/or 
biochemical characteristics including its pharmacokinetic properties in various anatomical and 
physiological compartments of the organism in question. The relevant physiological compartments are 
characterized with respect to the compound's various pharmacokinetic properties that may vary with the 
mode of administering the compound. A combined system of coupled differential equations is then used 
to simulate a concentration-time curve of the active compound in selected compartments. 
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(54) Bezeichnung: Verfahren zur Simulation der Wechselwirkung von chemischen Substanzen mit lebenden Orga- 
nismen 

(57) Zusammenfassung: Es wird ein Verfahren zur Ermitt- 
lung der Wechselwirkung einer oder mehrerer chemischer 
wirksamer Substanzen mit Organismen beschrieben. Das 
Verfahren umfasst die Schritte: 

A) Auswahl mindestens einer Substanz und Analyse ihrer 
physiologischen und/oder biochemischen Eigenschaften, 
insbesondere aus der Reihe Lipophilie, Loslichkeit, Prote- 
inbindung, Molekulgrofie (ausgedruckt als Molekulgewicht 
oder-volumen), pKa-Wertim Falle von Sauren oderBasen, 
metabolische Degradierungsrate und kinetische Konstan- 
ten von aktiven Transportern, 

B) Auswahl anatomischer und physiologischer Komparti- 
mente des jeweiligen Organismus, ausgewahlt wenigstens 
aus: a) im Falle von Saugetieren: Lunge, periphere Organe, 
Blutwege, bevorzugt arterielle und/oder venose Blutwege, 
oder Blutfliissigkeit; b) im Falle von Pflanzen Xylem- und 
Phloemstrom und Wurzel und/oder Blatt und/oder Stangel, 

C) Beschreibung der Kompartimente in einem Rechenpro- 
gramm im Hinblick auf die Transportprozesse der Substanz 
von und zu den Kompartimenten, Abbauprozesse und 
Wirkprozesse durch Massentransportgleichungen und ki- 
netische Reaktionsgleichungen zur Beschreibung des je- 
weiligen Prozesses, 

D) Beschreibung der Transportprozesse und gegebenen- 
falls der Abbauprozesse der wirksamen Substanz im Orga- 
nismus, insbesondere bei intravenoser, peroraler oder sub- 
kutaner Verabreichung, in dem Rechenprogramm durch 
Massentransportgleichungen und kinetische Reaktions- 
gleichungen, 

E) Verkniipfung der nach Schritt C) und □) ... 
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Beschreibung 

[0001] Die Erfindung betrifft ein Computerverfahren zur Berechnung des pharmakokinetischen und pharma- 
kodynamischen Verhaltens chemischer Substanzen in lebenden Organismen, z.B. in Saugetieren, Insekten 
oder Pflanzen. Die Raten und Ausmafie von Aufnahme, Verteilung, Metabolismus und Ausscheidung chemi- 
scher Substanzen in lebenden Organismen sind von groBer Bedeutung z.B. im Bezug auf die in vivo Aktivitat 
von pharmazeutischen Wirkstoffen und Pflanzenschutzmitteln oder fur die toxikologische Risikoabschatzung. 
Computerverfahren wie das physiologischbasierte pharmakokinetische (PB-PK) und pharmakodynamische 
(PD) Modelling sind geeignet, Wechselwirkungsprozesse von chemischen Substanzen mit lebenden Organis- 
men zu beschreiben. Die Erfindung ermoglicht eine besonders einfache Generierung solcher Modelle auf der 
Basis einer Multi-Ebenen Architektur, in der einzeln eingekapselte Module, die z.B. den Organismus, die Sub- 
stanz, die Applikationsform oder die Wirkform beschreiben, dynamisch miteinander verbunden werden. Durch 
Verbindung dieser Module zu einem Gesamtmodell wird ein System von gekoppelten Differentialgleichungen 
(Massenerhaltungsgleichungen) erzeugt, welches von einem Solver numerisch gelost wird. Der besondere 
Vorteil der Erfindung liegt in der hohen Flexibility des modularen Computerverfahrens, durch welches sich eine 
Vielzahl komplexer physiologischer und biochemischer Szenarien auf besonders einfache Weise generieren 
lasst. Die zur Erzeugung dieses Gleichungssystems erforderlichen chemischen, biologischen, physiologischen 
und anatomischen Input-Parameter konnen entweder in einer oder mehreren angebundenen Datenbanken 
enthalten sein oder iiber eine graphische Benutzeroberflache eingegeben werden. 

Stand der Technik 

[0002] Zahlreiche PB-PK- und PD-Modelle fur verschiedene Organismen wie Saugetiere [R. Kawai, M. Le- 
maire, J. L. Steimer, A. Bruelisauer, W. Niederberger, M. Rowland: Physiologically based pharmacokinetic stu- 
dy on a Cyclosporin derivate, SDZ IMM 125. J. Pharmacokin. Biopharm. 22, 327-365 (1994); P. Poulin, F.P. 
Theil: Prediction of pharmacokinetics prior to in vivo studies. 1 . Mechanism-based prediction of volume of dis- 
tribution. J. Pharm. Sci. 291, 129-156, (2002); P. Poulin, F.P. Theil: Prediction of Pharmacokinetics prior to in 
vivo studies. II. Generic physiologically based pharmacokinetic models of drug disposition, J. Pharm. Sci. 91, 
1358-1370 (2002)], Insekten [R. Greenwood, M. G. Ford, E. A. Peace, D. W. Salt: The Kinetics of Insecticide 
Action. Part IV: The in vivo Distribution of Pyrethroid Insecticides during Insect Poisoning. Pestic. Sci. 30, 
97-121 (1990)] oder Pflanzen [N. M. Satchlvi, E. W. Stoller, L. M. Wax, D. P. Briskin: A nonlinear dynamic si- 
mulation model for xenobiotik transport and whole plant allocation following foliar application. Pestic. Biochem. 
Physiol. 68, 67-84 (2000); N. M. Satchivl, E. W. Stoller, L. M. Wax, D. P. Briskin: A nonlinear dynamic simulation 
model for xenobiotik transport and whole plant allocation following foliar application. Pestic. Biochem. Physiol. 
68, 67-84 (2000)] sind in der Literatur beschrieben worden. Solche Modelle werden erfolgreich eingesetzt, urn 
das Verhalten chemischer Substanzen auf der Basis von physiologischen und anatomischen Informationen so- 
wie substanzspezifischen physiko-chemischen und biochemischen Parametern zu beschreiben. Die gangige 
Vorgehensweise dabei ist, den zu beschreibenden Organismus in mehrere physiologische Kompartimente zu 
unterteilen. Haufig entsprechen diesen Kompartimenten einzelne Organe wie Leber, Muskel, oder Lunge. Die- 
se Kompartimente sind durch physiologische Parameter wie z.B. Blutflussraten sowie Wasser-, Fett- und Pro- 
teinanteil charakterisiert [R. Greenwood, M. G. Ford, E. A. Peace, D. W. Salt: The Kinetics of Insecticide Action. 
Part IV: The in vivo Distribution of Pyrethroid Insecticides during Insect Poisoning. Pestic. Sci. 30, 97-121 
(1 990)]. Urn ein physiologie-basiertes Modell zu erhalten werden diese Kompartimente entsprechend den ana- 
tomischen Vorgaben des Organismus rechnerisch durch Massenerhaltungsgleichungen miteinander verbun- 
den, wie dies im oben zitierten Stand der Technik detailliert erlautert wird. 

Aufgabenstellung 

[0003] Nach dem Stand der Technik ist diese mathematische Kopplung starr, d. h. die Differentialgleichungen 
sind hart miteinander verdrahtet. Dadurch ist ein einmal erzeugtes Modell hinsichtlich der Zahl seiner Kompar- 
timente und deren physiologischer Verschaltung fixiert. Der Erfindung liegt die Aufgabe zu Grunde, einflexibles 
Verfahren zu entwickeln, welches eine variable Anzahl von Kompartimenten und Modulen (z.B. den Organis- 
mus, die Substanz oder die Applikationsform betreffend) nach individuellen Vorgaben dynamisch miteinander 
verschaltet, das entsprechende physiologisch-basierte Modell automatisch generiert und anschliefiend nume- 
risch lost. 

[0004] Gelost wird diese technische Aufgabe durch ein neuartiges modulares Programmierverfahren. Das 
Verfahren basiert auf einer Bibliothek aus voneinander unabhangigen Modulen, die jeweils fur sich abge- 
schlossen sind und verschiedene Funktionalitaten aufweisen (z.B. die Definition von Substanzeigenschaften 
oder Applikationsform oder die Beschreibung der Physiologie des zu untersuchenden Organismus, etc.). Erst 
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wahrend der Laufzeit des Programms werden diese Module dynamisch zu einem vollstandigen Modell zusam- 
mengestellt. Das hier beschriebene Verfahren erlaubt damit erstmalig eine flexible Behandlung unterschiedli- 
cher komplexer Szenarien ohne die Notwendigkeit, Anderungen auf der Ebene der Differentialgleichungen vor- 
zunehmen. 

[0005] Modulare Programmier-Konzepte sind vom Prinzip her aus anderen Anwendungsgebieten bekannt 
(siehe z.B. US 5,930,154 oderWO 00/65523), sie wurden aberbislang nichtfur pharmakokinetische oder phar- 
makodynamische Fragestellungen herangezogen. 

[0006] Gegenstand der Erfindung, durch die die genannte Aufgabe gelost wird, ist ein Verfahren zur Ermitt- 
lung der Wechselwirkung einer oder mehrerer chemischer wirksamen Substanzen mit Organismen mit den 
Schritten: 

A) Auswahl mindestens einer Substanz und Analyse ihrer physiologischen und/oder biochemischen Eigen- 
schaften, insbesondere aus der Reihe Lipophilie, Loslichkeit, Proteinbindung, MolekulgroUe (ausgedruckt 
als Molekulgewicht oder -volumen), pKa-Wert im Falle von Sauren oder Basen, metabolische Degradie- 
rungsrate und kinetische Konstanten von aktiven Transportern, 

B) Auswahl anatomischer und physiologischer Kompartimente des jeweiligen Organismus ausgewahlt we- 
nigstens aus: a) im Falle von Saugetieren: Lunge periphere Organe, Blutwege, bevorzugt arterielle 
und/oder venose Blutwege, oder Blutflussigkeit; b) im Falle von Pflanzen Xylem- und Phloemstrom und 
Wurzel und/oder Blatt und/oder Stangel 

C) Beschreibung der Kompartimente in einem Rechenprogramm im Hinblick auf die Transportprozesse der 
Substanz von und zu den Kompartimenten, Abbauprozesse und Wirkprozesse durch Massentransportglei- 
chungen und kinetische Reaktionsgleichungen zur Beschreibung des jeweiligen Prozesses 

D) Beschreibung der Transportprozesse und gegebenenfalls der Abbauprozesse der wirksamen Substanz 
im Organismus, insbesondere bei intravenoser, peroraler oder subkutaner Verabreichung, in dem Rechen- 
programm durch Massentransportgleichungen und kinetische Reaktionsgleichungen, 

E) Verknupfung der nach Schritt C) und D) ausgewahlten Gleichungssysteme und insbesondere numeri- 
sche Berechnung des entstehenden Systems gekoppelter Differentialgleichungen, 

F) Bestimmung des Konzentrations-/Zeitverlaufs der Wirksubstanz in ausgewahlten Kompartimenten. 

[0007] Bevorzugt steht der Organismus bzw. die Organismen Sauretiere, aus der Reihe Mensch, Affe, Hund, 
Schwein, Ratte oder Maus, oder Insekten, insbesondere Raupen, oder Pflanzen dar. 

[0008] Der Wirkstoff wird gernaR einer bevorzugten Ausfuhrung intravenos, oral, topikal, subkutan, intraperi- 
toneal, inhalativ uber die Nase oder Lunge, oder intracerebrovaskular (im Fall von Saugetieren), bzw. in- 
fra-haemolymphatisch, topikal, oder oral durch Fraft oder Schlundsonde (im Fall von Insekten) appliziert, oder 
die Applikation erfolgt als Spritzbruhe mit Aufnahme iiber Blatter oder Wurzeln (im Fall von Pflanzen). 

[0009] Bevorzugt ist ein Verfahren, das dadurch gekennzeichnet, dass die physiologischen Parameter, die die 
Organismen beschreiben, zeitabhangige Parameter sind. 

[001 0] In einer bevorzugten Ausfuhrung des Verfahrens werden die Parameter teilweise einer Datenbank ent- 
nommen werden, die an das Rechner-System angebunden ist. 

[0011] Die physiologischen und anatomischen Parameter konnen in einer bevorzugten Ausfuhrung mittels 
Zufallszahlen im Rahmen einer vorgegebenen statistischen Streuung variiert werden (Populationskinetik). 

[0012] Besonders bevorzugt ist ein Verfahren, bei denen fur das Rechenprogramm eine hierarchische Mo- 
dell-Architektur bestehend aus Modulen, die unterschiedliche Funktionalitaten aufweisen, und mindestens 
zwei Module, die entweder chemische oder biologische Eigenschaften der Substanz beinhalten, oder physio- 
logische oder anatomische Informationen uber den Organismus beinhalten, oder Informationen uber die Dar- 
reichungsform oder die Wirkart beinhalten, verwendet werden. 

[0013] Eine besonders bevorzugte Variante des Verfahrens ist dadurch gekennzeichnet, dass mehrere Mo- 
delle von Organismen zu einem Gesamtmodell, insbesondere zur Simulation von Drug-Drug-Wechselwirkun- 
gen, Mutter-Foetus-Modell, kombiniertes Insekt-Pflanzen-Modell, zusammengefasst werden. 

[0014] Kern dieser Erfindung ist die Realisierung eines modularen Simulationskonzepts mit dynamischer Ge- 
nerierung eines Differentialgleichungs- (DGL) Systems und dessen Anwendung auf pharmakokinetische und 
pharmakodynamische Fragestellungen. Das Simulationskernel besteht dazu aus einer Bibliothek einzelner 
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Module, die diefolgenden Funktionalitaten beinhalten: 

• Definition der physiko-chemischen Substanzeigenschaften (z.B. Lipophilie, Affinitat zu Plasmaproteinen, 
Molekulgewicht) 

• Anatomische und physiologische Beschreibung des Organismus (Saugetier, Insekt, Pflanze) 

• Beschreibung der einzelnen Kompartimente (Organe), die den Organismus bilden 

• Beschreibung der Applikationsform der chemischen Substanz (z.B. intravenos, oral, subkutan, topikal 
oder inhalativ fur Saugetiere, oral oder topikal fur Insekten, Applikation als Spritzbruhe und Aufnahme uber 
Blatter oder Wurzel fur Pflanzen) 

• Zeitschleifen-Definition zur numerischen Integration 

• Numerische Integration des Differentialgleichungssystems (z.B. nach dem Euler-Verfahren oder Run- 
ge-Kutta Verfahren) 

• Speichern der Simulationsergebnisse (z.B. in einer Datei, Datenbank, o. a.) 

[0015] Diese individuellen Module werden dynamisch wahrend der Laufzeit des Programms durch ein Com- 
putersystem, z.B. einen handelsublichen PC, auf eine vordefinierte Weise hierarchisch miteinander verbun- 
den. AnschlielJend wird automatisch ein System von gekoppelten DGLs (Massenerhaltungsgleichungen) er- 
zeugt, welches durch numerische Integration gelost wird. Als Ergebnis liefert das Modell Konzentrati- 
ons/Zeit-Kurven der betrachteten chemischen Substanz in den diversen Kompartimenten des betrachteten Or- 
ganismus, welche dann weiter z.B. zur Berechnung einer pharmakologischen Wirkung herangezogen werden 
konnen. 

[0016] Das minimale Modell, welches der Erfindung zu Grunde liegt, besteht mindestens aus einem Sub- 
stanz-Modul, einem Applikations-Modul, einem Organismus-Modul, dem DGL-Generator und dem numeri- 
schen Solver fur das DGL-System, sowie optionalen weiteren Modulen, z.B. einem pharmakodynamischen 
Wirkmodul. 

[0017] Das Substanz-Modul beinhaltet physiko-chemische und/oder biochemische Informationen uber die 
Substanz, deren Verhalten simuliert werden soil. Dies konnen z.B. Werte fur die Lipophilie, Loslichkeit, z.B. 
Wasser oder intestinalem Fluid, Proteinbindung, z.B. von Plasmaproteinen, Molekulgrofie (ausgedruckt als 
Molekulgewicht oder -volumen), pKa-Wert im Falle von Sauren oder Basen, metabolische Degradierungsra- 
ten, kinetische Konstanten von aktiven Transportern, etc. sein. Diese Parameter sind generell durch in vitro 
Experimente bestimmbar oder - in einzelnen Fallen, z.B. im Fall der Lipophilie - durch bekannte Vorhersage- 
modelle, die z.B. QSAR, HQSAR oder neuronalen Netze verwenden, direkt aus der Struktur berechenbar. 

[0018] Das Organismus-Modul beinhaltet physiologische und anatomische Informationen, die den Organis- 
mus charakterisieren, mit dem die Substanz wechselwirken soli. Organismus-spezifische Informationen sind 
unter anderem Volumina und Volumenblutflussraten sowie Wasser-, Fett- und Proteingehalt der einzelnen 
Kompartimente (im Fall von Saugetieren und Insekten), oder Xylem- und Phloemtransportraten und die Volu- 
mina von Symplast, Apoplast und Vakuolen im Fall von Pflanzen. Diese physiologischen Parameter konnen 
zeitlich konstant sein oder auch zeitabhangig sein, urn z.B. Wachstumsprozesse und andere physiologische 
Veranderungen mit der Zeit mit berucksichtigen zu konnen. Diese physiologischen und anatomischen Para- 
meter sind fur eine Vielzahl relevanter Organismen in der Literatur veroffentlicht. 

[0019] Das Applikations-Modul beinhaltet Informationen uber den Ort und das zeitliche Profil der Verabrei- 
chung der Substanz. Ubliche Formen fur die Verabreichung von pharmazeutischen Wirkstoffen bei Saugetie- 
ren sind die intravenose Applikation als Bolus oder Infusion, orale Applikation als Losung, Kapsel oder Tablette, 
subkutane Administration, intraperitoneale Administration, topikale Administration mit Aufnahme uber die Haut 
oderSchleimhaute, sowie nasaleoder inhalative Administration. Insekten nehmen chemische Substanzen vor- 
wiegend oral oder uber die Kutikula nach topikalem Kontakt auf. Die Substanzaufnahme in Pflanzen verlauft 
entweder uber die Wurzeln oder uber die Blatter. 

[0020] Der DGL-Generator erzeugt automatisch des Differentialgleichungssystem basierend auf der Massen- 
bilanz innerhalb des Organismus. Der Solver ubernimmt die numerische Integration des DGL-Systems und lie- 
fert als Resultat Konzentrations-Zeit-Profile der chemischen Substanz in den Kompartimenten des Organis- 
mus. Diese Daten konnen weiter genutzt werden urn z.B. eine pharmakologische Wirkung an einem Zielenzym 
(Target) zu beschreiben. 

[0021] Zur dynamischen Verknupfung der Module ist ein hierarchische Management der Variablen, die in den 
Modulen enthalten sind, erforderlich, da verschiedene Variablen aus unterschiedlichen Modulen voneinander 
abhangen konnen. So werden z.B. die substanz-spezifischen physiko-chemischen Parameter (aus dem Sub- 
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stanz-Modul) zusammen mit physiologischen Informationen (aus dem Organismus-Modul) verwendet, um 
Gleichgewichts-Verteilungskoeffizienten zwischen Plasma und den peripheren Kompartimenten im Saugetier- 
korperzu berechnen. In ahnlicher Weise existieren Abhangigkeiten zwischen den einzelnen Kompartimenten 
innerhalb eines Organismus. Der Blutfluss in der Lunge ist beispielsweise gegeben durch die Summe der Blut- 
flussraten aller ubrigen Organe, da die Blutzirkulation im Saugetierkorper in sich geschlossen ist. Solche Ab- 
hangigkeiten erfordern eine generelle hierarchische Datenstruktur, damitsie automatisch erkannt und beruck- 
sichtigt werden konnen. Eine solche hierarchische Datenstruktur ist Teil der vorliegenden Erfindung. Die Mo- 
dule werden hierzu von einer Datenbank von Objekten verwaltet. Die Objektdatenbank kann dynamisch belie- 
big erweitert werden, was ein Hochstmali an Funktionalitat garantiert. Auch mehrere Module mit gleicher oder 
ahnlicher Funktionalitat sind zulassig. Die Datenstrukturen fur die mathematische Beschreibung des Gesamt- 
modells sind in den einzelnen Modulen eingekapselt. Eine Besonderheit dieses Konzeptes ist die moduluber- 
greifende hierarchische Verknupfung der Daten, die wiederum, ahnlich den Objektmodulen, von einer eigenen 
Datenbank verwaltet werden. Dadurch ist ein, modellbedingter, Zugriffauf Daten anderer Module moglich. Die 
Module des Simulationskernels mussen von einem Hauptprogramm generiert und in die gewunschte Reihen- 
folge verknupft werden. 

[0022] Zur Verwaltung der diversen Input-Parameter wurde zum einen eine Datenbankanbindung der Com- 
pound- und Organismus-Parameter realisiert. Daruber hinaus (oder auch als einzige Alternative) stellt ein gra- 
phisches User-Interface (GUI) sicher, dass alle relevanten Model l-Parameter (a) fur den Usersichtbar und (b) 
durch ihn editierbarsind. Das Verandern einzelner Daten in einem Modul fuhrt zur Prufung der Abhangigkeiten 
in den restlichen, entsprechend der Datenhierarchie, uber vordefinierte Meldungen. Diese Prufung der Da- 
ten-Abhangigkeiten ist wesentlicher Bestandteil des Modell-Generators, da sie zu jedem Zeitpunkt die Korrekt- 
heit des Modells sicherstellt. Uber das GUI werden auch die Ergebniskurven sowie optional aus diesen Kurven 
abgeleitete typische pharmakokinetische Parameter (z.B. Flache unter der Kurve, maximale Konzentration, 
Zeitpunkt der maximalen Konzentration, Halbwertszeit, etc.) oder pharmakodynamische Parameter (z.B. Wirk- 
starke und -dauer) angezeigt. 

[0023] Der besondere Vorteil der vorliegenden Erfindung liegt in der Flexibility, mit der auf einfache Weise 
unterschiedliche und komplexe Modelle erzeugt werden konnen. Diefolgenden Beispiele sollen dies anhand 
von Figuren illustrieren. Sie zeigen relevante Szenarien unterschiedlicher Komplexitat. Die Beispiele sollen 
nicht als Limitierung der Anwendbarkeit der vorliegenden Erfindung verstanden werden. 

Ausfuhrungsbeispiel 



[0024] 


In den Figuren zeigen: 






[0025] 


Fig. 


1 Die schematische 


Darstellung 


eines vier-Kompartiment Saugetiermodells 


[0026] 


Fig. 


2 Die schematische 


Darstellung 


eines physiologischen Ganzkorper-Saugetiermodells 


[0027] 


Fig. 


3 Die schematische 


Darstellung 


eines Organs im Ganzkorpermodell 


[0028] 


Fig. 


4 Die schematische 


Darstellung 


eines physiologischen Modells fur Raupen 


[0029] 


Fig. 


5 Die schematische 


Darstellung 


eines physiologischen Pflanzenmodells 


[0030] 


Fig. 


6 Die schematische 


Darstellung 


des einfachsten Modell-Szenarios: Das Modell besteht aus den 



Modulen "Substanz", "Verabreichung", "Organismus", "DGL-Generator", und "Integrator". Optional kann ein 
"Aktions-Modul" betrachtet werden. 

[0031] Fig. 7 Die schematische Darstellung der Anbindung des Simulationsmodells an Datenbanken und die 
graphische Benutzeroberflache. 

[0032] Fig. 8 Die schematische Darstellung eines Modells fur multiple Administration 

[0033] Fig. 9 Die schematische Darstellung eines Modells fur eine Substanz, die einem enterohepatischen 
Kreislauf unterliegt 

[0034] Fig. 10 Die schematische Darstellung eines Modells zur Beschreibung der Wechselwirkung zweier 
Substanzen 
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[0035] Fig. 11 Die schematische Darstellung eines Modells zweier miteinander in Wechselwirkung stehender 
Organismen. Beispiele fur ein solches Szenario sind des Mutter-Foetus-Modell sowie das kombinierte Insek- 
ten-Pflanze-Modell. 

[0036] Fig. 12 den Konzentrations/Zeitverlauf einer Substanz in verschiedenen Organen einer Ratte 

[0037] Fig. 13 den Konzentrations/Zeitverlauf einer Substanz im Plasma verschiedener Sauger 

[0038] Fig. 14 den Konzentrations/Zeitverlauf einer Substanz in verschiedenen Organen eines Menschen 

[0039] Fig. 15 den Konzentrations/Zeitverlauf einer Substanz nach peroraler Administration im Plasma eines 
Menschen 

Beispiele 

[0040] In diesem Abschnitt werden bevorzugte Ausfuhrungsformen fur die zu beschreibenden Organismen 
(Saugetier, Insekt und Pflanze) beschrieben, anschlieftend wird deren Einbau in das dynamische Simulations- 
modell an Beispielen gezeigt. 

Beispiele fur verschiedene Organismus-Module 

[0041] Fig. 1 zeigt ein stark vereinfachtes vier-Kompartiment-Modell fur einen Saugetierkorper. Der Sauge- 
tierkorper besteht hier lediglich aus venosen und arteriellen Blutpool, der Lunge und einem weiteren periphe- 
ren Kompartiment, uberdas die Substanz metabolisch eliminiert wird. Das periphere Kompartiment wird iiber 
den arteriellen Blutstrom (Blutflussrate Q org ) mit der Substanz (Blutkonzentration C in org ) versorgt. Nach Wech- 
selwirkung in dem Kompartiment stromt das venose Blut (Konzentration der Substanz C out org ) in die Lunge, iiber 
die der Blutkreislauf geschlossen ist (Q| Unge = Q org ). 

[0042] In einer bevorzugten Ausfuhrungsform der Erfindung wird als Organismus-Modul das in Fig. 2 gezeig- 
te Saugetiermodell verwendet. Der Saugetierkorper besteht hier ebenfalls aus dem arteriellen und dem veno- 
sen Blutpool, einem Lungen-Kompartiment und den folgenden peripheren Organen: Leber, Niere, Muskel, 
Knochen, Haut, Fett, Gehirn, Magen, Dunndarm, Dickdarm, Bauchspeicheldruse, Milz, Gallenblase und Ho- 
den. Diese Ausfuhrungsform wird im folgenden als "Ganzkorper-Modell" bezeichnet. Jedes Organ In diesem 
Ganzkorper-Modell ist wiederum unterteilt in mehrere Subkompartimente, die jeweils den vaskularen Raum 
(bestehend aus Plasma und den roten Blutkorperchen), interstitiellen und intrazellularen Raum reprasentieren 
(Fig. 3). Plasma und interstitieller Raum werden als equilibriert behandelt. Der Substanztransport zwischen 
dem interstitiellen Raum und dem intrazellularen Raum wird beschrieben als passiver Diffusionsprozess, der 
einer Kinetik erster Ordnung folgt, oder als aktiver Transportprozess, der mittels einer sattigbaren Michae- 
lis-Menten Kinetik modelliert wird. Analog existieren in jedem intrazellularen Subkomprtiment ein oder mehrere 
Michaelis-Menten-Terme, die dem metabolischen Abbau der Substanz Rechnung tragen. 

[0043] Das resultierende System von gekoppelten Differentialgleichungen hat die folgende Form: Fur jedes 
periphere Organ (oberer Index "org") existieren drei Differentialgleichungen fur die Subkompartimente "Plas- 
ma" (pi), "Rote Blutkorperchen" bzw. "Blutzellen" (be) und "Zellinneres" (cell). Das Ganzkorpermodell beinhal- 
tet die folgenden peripheren Organe: Lunge, Magen, Dunndarm, Dickdarm, Pankreas, Milz, Leber, Niere, Ge- 
hirn, Herz, Muskel, Knochen, Haut, Fett und Hoden. Fur das Plasma und den interstitiellen Raum (int) ergibt 
sich pro Organ: 
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MASSENBILANZGLEICHUNG BEDEUTUNG 

dC°, s Inter-kompartimenteller 

[C(l-HCT) + C]V OTg -^ 

dt =Q7(qr-c?) 



Diffusiver 

Massentransport zu den 
roten Blutkorperchen 
Diffusiver 



PA org I P org 

_ r ^ porg _ ^cell 

K. , I K org J Massentransport 

Zellirmere 

V org C°[ s / K , Aktiver Transportterm in 

j£ org . _)_ / K , ^ as Zellirmere hinein 

(Influx) 

yorg Qoig j ^jT-org jr- ^ ^ Aktiver Transportterm aus 
+ K° rg + C° r \ I (K° rs K j) ^em Zellinneren heraus 
(Efflux) 

[0044] Fur die roten Blutkorperchen folgt: 

MASSENBILANZGLEICHUNG BEDEUTUNG 

An ot % Inter-kompartimenteller 

Cf HCT V OT * 

fa _ Qorg ^qzti _ Flussterm 



Diffusiver 

Massentransport zu i 
roten Blutkorperchen 



[0045] Das Zellinnere wird beschrieben mittels: 
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MASSENBILANZGLEICHUNG 



riC crg 

dt 



PA org ( porg ^ 

. rA | r org _ '-"cell 



VZ,exCW(K org K pl ) 
K^x +C° e r f,/(K ors K ,) 



2005.04.21 

BEDEUTUNG 

Diffusiver Massen- 
transport 



max,Ml cell 



CS/(K-»K„) 



org porg / /-tj-org \ 

m.Ml + '-cell ' ^p\) 

T Z,M2 C^/(K^K pl ) 



K£ M2 + CS / (K or 

CL org porg 



K P.) 



Zellinnere 

Aktiver Transportterm 
in das Zellinnere 
hinein (Influx)) 
Aktiver Transportterm 
aus dem Zellinneren 
heraus (Efflux) 
Metabolischer Abbau 
(Enzym 1) 

Metabolischer Abbau 
(Enzym 2) 

Abbau erster Ordnung 
(z.B. in der Leber oder 
Niere) 



[0046] Die Lunge (org = Ing) wird analog beschrieben , allerdings ist hier der Zustrom nicht uber den arteriellen 
Blutpool (art) sondern uber den venosen (ven): 



MASSENBILANZGLEICHUNG 

dC lng 

[C s g (l-HCT) + C g ]V lng - 



BEDEUTUNG 

Blurfluss ist umgekehrt im 
Vergleich zu den iibrigen 
Organen ! 

analog den anderen 
Organen 

Blutfluss ist umgekehrt im 
Vergleich zu den iibrigen 
Organen ! 

± ... analog den anderen 

Organen 

[0047] Die aus den Organen herausstromenden Konzentrationen vereinigen sich im venosen Blutpool. Die 
resultierende Plasmakonzentration folgt aus dem blutfluss-gewichteten Mittelwert der einzelnen Organkonzen- 
trationen: 



f '" s HCT V lng 
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MASSENBILANZGLEICHUNG 

dt 



BEDEUTUNG 



Diffusiver 
Massentransport 



dt 

x _ J 1 - #Cr /or jotewza 
~ |#CT /or red Wood cefts 

[0048] Der arterielle Blutpool wird schlieRlich beschrieben durch: 

MASSENBILANZGLEICHUNG 



Blutkorperchen (+) 
Plasma Clearance (ggf.) 

Input-Funktion fur 
Application in den 
venosen Blutpool 



BEDEUTUNG 



lv "pl/bc 
dt 



= Qpl/bc (Cpl/1 

_ pa; 



Cpl/bc) 



Diffusiver Massentrans- 
port zwischen Plasma 
(-) und den roten Blut- 
korperchen (+) 




[ l-HCT for plasma 
HCT for red blood cells 



[0049] In diesem Gleichungssystem bedeuten: 

f x ° rg : Volumenanteil des vaskularen (x = vas), interstitiellen (x = int) oder intrazellularen Raums (x = cell) des 
Organs 

V° r9 : Volumen des Organs 

HCT: Haematokrit (= Volumenanteil der roten Blutkorperchen im Vollblut) 
Q x ° rg : Blutflussrate des Plasmas (x = pi) bzw. der Blutzellen (x = be) im Organ 
PA bc ° rg : Permeabilitats-Oberflachenprodukt der roten Blutkorperchen 
PA °rg. Permeabilitats-Oberflachenproduktdes Organs 

K p ,: Gleichgewichts-Verteilungskoeffizient der Substanz zwischen Plasma und Wasser 
K bc : Gleichgewichts-Verteilungskoeffizient der Substanz zwischen roten Blutkorperchen und Plasma 
K° rg : Gleichgewichts-Verteilungskoeffizient der Substanz zwischen dem Organ und Plasma 
V max , K m : Michaelis-Menten-Konstanten fur aktiven Transport bzw. Metabolismus 



[0050] Alternativ zu oder in Kombination mit der hier beschriebenen intravenosen Applikation kann auch eine 
Aufnahme uberdie Magen-Darm-Schleimhaut nach peroraler Verabreichung simuliert werden. Die Losungen 
dieses Gleichungssystems ergeben die Konzentrations-Zeit-Beziehungen furalle im Modell enthaltenen Kom- 
partimente. 
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[0051] Zur Beschreibung einer pharmakologischen Wirkung kann weiterhin die Konzentrations-Zeit-Bezie- 
hung in dem Kompartiment, welches das biologische Target des Wirkstoffs enthalt, mit einem pharmakodyna- 
mischen Effekt verknupft werden. Typische Effektfunktionen sind z.B.: 

• Hyperbolische oder sigmoide Emax-Modelle: 

Effekt = E 0 + Emax C * 
EC} 0 + C{ 

Effekt = pharmakologischer Wirkparameter (zeitabhangig) 

E 0 = Basiswert des pharmakologischen Wirkparameters 

E max = Maximalwert der pharmakologischen Wirkung 

EC 60 = Konzentration, bei der 50 % des maximalen Effekts erreicht sind 

C x = Konzentration am Wirkort (zeitabhangig) 

Y = Formparameter 

• Potenzfunktionen: Effekt = E 0 + |3 C x bzw. Log-Linear Modelle: 
Effekt = E 0 + p Ln(C x ) 

Effekt = pharmakologischer Wirkparameter (zeitabhangig) 

E 0 = Basiswert des pharmakologischen Wirkparameters 

P = Parameter fur die Steigung des Effekts als Funktion der Konzentration 

C x = Konzentration am Wirkort (zeitabhangig) 

Y = Formparameter 

• Wirkstoffinteraktionsmodelle wie z.B. partieller oder vollstandiger Antagonismus, etc. 

• Kombinationen der vorgenannten Modelle, mit denen z.B. multiple Wirkzentren, oder Rezeptor-Transdu- 
cer-Wechselwirkungen beschrieben werden konnen. 

[0052] In einer weiteren bevorzugten Ausfuhrungsform der Erfindung reprasentiert das Organismus-Modul 
den Korperbau eines Insekts, insbesondere den einer Raupe (Fig. 4). Der Raupenkorper besteht aus den fol- 
genden Kompartimenten: Hamolymphe als zentrales Kompartiment, Kutikula, Muskel, Fettkorper, Nervensys- 
tem, und Darmwand als periphere Kompartimente, sowie den Kompartimenten Kutikula-Oberflache und 
Darminhalt, Liber die ein Substanzaustausch mit der Umgebung stattfinden kann. Der interkompartimentelle 
Massentransport kann wiederum passiv via Diffusion oder mit Hilfe von Transportern als aktiver Prozess erfol- 
gen. Das resultierende Differentialgleichungssystem ist beschrieben worden im Anhang zu der noch nicht ver- 
offentlichten deutschen Patentanmeldung mit dem Aktenzeichen 10256315.2. 

[0053] Das Pflanzenmodell aus Fig. 5 stellt eine weitere bevorzugte Ausfuhrungsform dar. Gemafi der typi- 
schen Pflanzenphysiologie beschreibt das Modell die Wurzeln, Stangel, und Blatter einer Pflanze. Jedes dieser 
Kompartimente besteht seinerseits aus drei Subkompartimenten, welche die Vakuole, den Symplasten und 
den Apoplasten reprasentieren. Diese Subkompartimente sind untereinander durch Membranen voneinander 
separiert. Diese Membranen konnen - wie in den vorab beschriebenen Organismen - auch hier durch passive 
Diffusion oder mittels aktivem Transport permeiert werden. Weiterhin sind die Subkompartimente durch unter- 
schiedliche pH-Werte gekennzeichnet, welche insbesondere auf das Verteilungsverhalten von Sauren und Ba- 
sen einen starken Einfluss ausuben. Zwischen den Kompartimenten verfugt die Pflanze iiber zwei Transport- 
pfade: Der Xylemstrom im Apoplasten fliefit von den Wurzel in Richtung Blatter, wahrend der Phloemstrom um- 
gekehrt Stoffe im Symplasten von den Blattern zu den Wurzeln transportiert. 

Beispiele fur Modellstrukturen 

• Die einfachste Modellstruktur (Fig. 6) 

[0054] Fig. 6 zeigt die einfachste Modellstruktur. Das komplette Modell besteht aus den Modulen "Substanz", 
"Verabreichung", "Organismus", "DGL-Generator", und "Integrator". Optional kann ein "Aktions-Modul" be- 
trachtet werden. Fig. 7 zeigt schematisch die Anbindung an eine oder mehrere Datenbanken, die z.B. die Sub- 
stanzparameter oder die physiologischen und anatomischen Informationen des betrachteten Organismus be- 
inhalten konnen. 

• Multiple Administration (Fig. 8) 
[0055] Die wiederholte Verabreichung eines Wirkstoffs an ein und denselben Organismus ist ein realistisches 
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Szenario fur eine Vielzahl von pharmazeutischen Substanzen, die in regelmaftigen Abstanden uber einen lan- 
geren Zeitraum genommen werden mussen wie z.B. antiinfektive Wirkstoffe. Mit diesem Beispielmodell kon- 
nen Langzeiteffekte wie z.B. die Akkumulation des Wirkstoffs in einzelnen Organen simuliert werden. 

• Enterohepatischer Kreislauf (Fig. 9) 

[0056] Eine Reihe von chemischen Stoffen, die an Saugetiere verabreicht wurden, unterliegen einem entero- 
hepatischen Kreislauf. In diesem Fall wird ein Teil der verabreichten Substanz unverandert in der Leber in die 
Gallenflussigkeit abgesondert und darin in der Gallenblase (reprasentiert durch ORGAN N in Fig. 9) akkumu- 
liert. Ausgelost durch einen chemischen Reiz, der z.B. durch die Einnahme einer Mahlzeit getriggert werden 
kann, zieht sich die Gallenblase zusammen und gibt ihren Inhalt - Gallenflussigkeit, die unter anderem auch 
den Wirkstoff enthalt - in den Zwolffingerdarm ab. Uber den Darm kann die Substanz denn absorbiert werden 
und somit erneutTeil dersystemischen Zirkulation werden. Im in Fig. 8 skizzierten Modell kann dieser Prozess 
einfach als eine erneute intestinale Verabreichung (APPLICATION 2, charakterisiert durch eine Lag-Zeit) eines 
Teils der ursprunglichen Substanz beschrieben werden. 

• Drug-Drug Wechselwirkung (Fig. 10) 

[0057] Wechselwirkungen von mehreren Substanzen (in diesem Beispiel: zwei Substanzen) untereinander 
sind von grafter Bedeutung. Modellhaft kann dieser Fall wie folgt abgebildet werden: ORGANSIM 1 und 2 sind 
uber identische physiologische Parameter definiert, denn sie reprasentieren ein und denselben Organismus. 
Zwei unterschiedliche Substanzen werden separat verabreicht, die Art, Zeitpunkt, und Dauer der Verabrei- 
chung kann fur beide Substanzen unterschiedlich sein. Die eigentliche Wechselwirkung der beiden Substan- 
zen wird in einem Aktionsmodul definiert, welches die Organe ORGAN N in den Organismen 1 und 2 verknupft. 
Die Wechselwirkung kann beispielsweise eine kompetitive Inhibition oder jede andere denkbare biochemische 
Interaktion sein. 

• Mutter-Foetus Model (Fig. 11) 

[0058] Das Mutter-Foetus Modell in Fig. 11 ist ein Beispiel, bei dem zwei unterschiedliche aber gekoppelte 
Organismen simultan modelliert werden. ORGANISM 1 reprasentiert die Mutter und ORGANISM 2 den Foe- 
tus. Die Blutzirkulation des Feten ist an den Blutkreislauf der Mutter uber die Plazenta (hier reprasentiert durch 
ORGAN N in ORGANISM 1) angeschlossen. In diesem Beispiel ist von grafter Bedeutung, dass die physiolo- 
gischen Parameter, insbesondere die Organvolumina und Blutflussraten des Feten, Funktionen der Zeit sein 
konnen, urn das genaue Gestationsalter sowie das Wachstum des Feten im Mutterleib berucksichtigen zu kon- 
nen. 

• Kombinierte Pflanzen/lnsekt Modelle (Fig. 11 ) 

[0059] Ahnlich wie das Mutter-Foetus-Model I stellt auch das Pflanze/lnsekt Modell eine Kombination zweier 
gekoppelter Organismen dar, allerdings sind die Organismen hier physiologisch unterschiedlich. ORGANISM 
1 beschreibt eine Pflanze, die eine applizierte Substanz, z.B. ein Insektizid, uber die Blattkutikeln oder die Wur- 
zeln aufnimmt. Nach Verteilung im gesamten Organismus steht diese Substanz dann auch in den ubrigen 
Kompartimenten zur Verfugung. Durch Fraft an einem Blatt, welches ein spezielles Kompartiment in dem 
Pflanzen-Organismus darstellt, kann die Substanz von einem Insekt (ORGANISM 2) aufgenommen werden, 
verteilt sich dort und erzielt schlieftlich im Zielorgan, beispielsweise dem Nervensystem der Raupe, seine In- 
sektizide Wirkung. 

• Weitere Kombinationen dervorab beschriebenen Beispiele 

[0060] Kombinationen dervorab beschriebenen Beispiele sind ebenfalls von grafter Bedeutung. Die multiple 
Administration von zwei oder mehreren miteinander wechselwirkenden Substanzen in einem Mutter-Foe- 
tus-Modell z.B. kann zur fruhzeitigen Abschatzung des toxikologischen Risikos des Feten herangezogen wer- 
den. 

• Simulationsergebnisse fur das Ganzkorper-Saugetiermodell 

[0061] Im folgenden wird beispielhaft gezeigt, welche Schritte zur Durchfuhrung einer Simulation erforderlich 
sind. Zunachst mussen die substanzabhangigen Eigenschaften bestimmt werden. Als Beispiel dient hier eine 
Substanz X mit den folgenden Eigenschaften: 
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Parameter 


Wert 


Einheit 


Lipophilie (LogMA) 


2.7 


-/- 


ungebundene Plasmafraktion 


0.25 


-/- 


Leberclearance 


1 


ml/min/kg 


Dosis 


1 


mg/kg 



Tabelle 1 : Substanzabhangige Parameter der Substanz X 

[0062] Aus der Lipophilie (MA = 10 A LogMA) und der ungebundenen Plasmafraktion (f u ) konnen dann mit Hilfe 
von veroffentlichten Gleichungen die Organ-Verteilungskoeffizienten berechnet werden [M. Harter, J. Kelde- 
nich, W. Schmitt: Estimation of physicochemical and ARME parameters, Ch. 26 in: Combinatorial Chemisrty- 
A Practical Handbook, Part IV. Eds. K. C. Nicolau et al., Wiley VCH, Weinheim, 2002]. Fur die Organ-Vertei- 
lungskoeffizienten ergeben sich in diesem Beispiel die folgenden Werte: 



Organ/Plasma Verteilungskoeffizient 

______ __ 

Dunndarm 8.34 



Dickdarm 8.34 

Pankreas 10.57 

Milz 2.74 

Leber 9.35 

Niere 7.19 

Lunge 1.97 

Gehirn 14.21 

Herz 13.28 

Muskel 2.33 

Knochen 34.45 

Haut 13.49 

Fett 100.42 

Hoden 4.42 



Tabelle 2: Aus den Werten von Tabelle 1 resultierende Verteilungskoeffizienten [8]. 

[0063] Weiterhin mussen die physiologischen Parameter wie Organvolumina, -blutflussraten und Zusammen- 
setzung hinsichtlich vaskularem, interstitiellen und zellularem Raum bekanntsein. Diese sind ebenfalls in der 
Literatur beschrieben. Fur die Spezies Maus, Ratte, Hund und Mensch finden sich beispielsweise die in Tabel- 
len 3 bis 5 aufgelisteten Werte: 
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Organvolumina [ml] 


MAUS 


RATTE 


HUND 


MENSCH 


venoser Blutpool 


0.128 


6.8 


0.8 


250 


arterieller Blutpool 


0.073 


2.9 


0.8 


140 


Lunge 


0.201 


2.2 


145 


670 


Magen 


0.1 


1.1 


40 


150 


Diinndarm 


0.2 


11.1 


140 


640 


Dickdarm 


0.2 


11.1 


140 


370 


Pankreas 


0.13 


1.3 


10 


100 


MHz 


0.13 


1.3 


10 


180 


Leber 


0.941 


10 


366 


1710 


Gallenblase 


0.12 


1.2 


10 


20 


Niere 


0.449 


7 


154 


720 


Gehirn 


0.336 


1.671 


90.4 


1486 


Herz 


0.5 


1.2 


90 


330 


Muskel 


11.74 


110.1 


6502 


30200 


Knochen 


2.775 


28.2 


2584 


12060 


Haut 


5.16 


43.4 


647 


3020 


Fett 


1.51 


14.2 


2670 


10060 


Hoden 


0.001 


2.5 


5 


35 


Summe der Organe 


24.694 


257.271 


13605 


62141 



Tabelle 3: Organvolumina fur Maus, Ratte, Hund und Mensch (Literaturdaten) 



Blutfluss [ml/min] 


MAUS 


RATTE 


HUND 


MENSCH 


Magen 


0.5 


1.22 


40 


60 


Diinndarm 


0.5 


7.02 


400 


600 


Dickdarm 


0.5 


3.82 


160 


240 


Pankreas 


0.005 


0.51 


40 


60 


Milz 


0.063 


0.63 


120 


180 


Leber 


2.25 


5.5 


310 


390 
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Blutfluss [ml/min] 


MAUS 


RATTE 


HUND 


MENSCH 


Niere 


0.67 


14.6 


243 


1133 


Gehirn 


0.11 


1.1 


145 


700 


Herz 


0.039 


3.92 


180 


240 


Muskel 


0.33 


7.2 


118 


550 


Knochen 


0.007 


1.6 


35 


167 


Haut 


0.02 


4.8 


10 


50 


Fett 


0.002 








Hoden 


0.0002 


0.48 


0.66 


2.6 


Lunge (=Summe) 


4.996 


53.72 


1805 


4670 



Tabelle 4: Organblutflussraten fur Maus, Ratte, Hund und Mensch (Literaturdaten) 



Volumenanteil 


f_ 


_vas 


f 


Jnt 


f_cell 


Fettgewebe 


0 


.010 


0 


.135 


0.855 


Gehirn 


0 


.037 


0 


.004 


0.959 


Magen-Darm-Trakt 


0 


.032 


0 


.100 


0.868 


Herz 


0 


.262 


0 


.100 


0.638 


Niere 


0 


,105 


0 


.200 


0.695 


Leber 


0 


.115 


0 


.163 


0.722 


Lunge 


0 


,626 


0 


188 


0.186 


Muskel 


0 


,026 


0, 


120 


0.854 


Knochen 


0 


041 


0 


100 


0.859 


Haut 


0. 


,019 


0. 


302 


0.679 


Bauchspeicheldriise 


0. 


180 


0. 


120 


0.700 


Milz 


0, 


282 


0, 


150 


0.568 


Hoden 


0, 


140 


0. 


069 


0.791 



Tabelle 5: Zusammensetzung der Organe von Saugetieren (Literaturdaten) 

[0064] Im Folgenden sind Simulationsergebnisse gezeigt, die mit dem Ganzkorpermodell nach Fig. 3 erzielt 
wurden. Hierbei wurde als einfachster Grenzfall angenommen, dass der Transport in die Organe blutflusslimi- 
tiertist (d. h. PA° rg ->«,). 

[0065] Fig. 12 zeigt die resultierenden Organkonzentrationen in der Ratte nach einer intravenosen Bolus Ad- 
ministration von 1 mg/kg der Substanz X. 

[0066] Fig. 13 zeigt die Konzentration-Zeit-Kurve im Plasma simuliert in den vier verschiedenen Spezies 
nach einer intravenosen Bolus Administration von 1 mg/kg der Substanz X. 

[0067] Fig. 14 zeigt die Konzentrations-Zeit Kurven in den diversen Organen eines Menschen nach mulitpler 
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intravenoser Bolus Administration von 1 mg/kg der Substanz X uber drei Tage (Dosierungsschema: 6 h - 6 h 
- 12 h, berechnet nach dem Modell in Fig. 8) 

[0068] Fig. 15 zeigt Konzentration-Zeit-Kurve im Plasma nach peroraler Administration von 1 mg/kg der Sub- 
stanz X im Menschen ohne enterohepatischem Kreislauf und mit enterohepatischem Kreislauf berechnet nach 
dem Modell in Fig. 9 (der Anteil der biliaren Clearance an dertotalen ist mit 25 % angenommen). 

Patentanspriiche 

1 . Verfahren zur Ermittlung der Wechselwirkung einer oder mehrerer chemischer wirksamen Substanzen 
mit Organismen mit den Schritten: 

G) Auswahl mindestens einer Substanz und Analyse ihrer physiologischen und/oder biochemischen Eigen- 
schaften, insbesondere aus der Reihe Lipophilie, Loslichkeit, Proteinbindung, MolekulgroUe (ausgedruckt als 
Molekulgewicht oder -volumen), pKa-Wert im Falle von Sauren oder Basen, metabolische Degradierungsrate 
und kinetische Konstanten von aktiven Transportern, 

H) Auswahl anatomischer und physiologischer Kompartimente des jeweiligen Organismus ausgewahlt wenigs- 
tens aus: a) im Falle von Saugetieren: Lunge periphere Organe, Blutwege, bevorzugt arterielle und/oder veno- 
se Blutwege, oder Blutflussigkeit; b) im Falle von Pflanzen Xyiem- und Phloemstrom und Wurzel und/oder Blatt 
und/oder Stangel 

I) Beschreibung der Kompartimente in einem Rechenprogramm im Hinblick auf die Transportprozesse der 
Substanz von und zu den Kompartimenten, Abbauprozesse und Wirkprozesse durch Massentransportglei- 
chungen und kinetische Reaktionsgleichungen zur Beschreibung des jeweiligen Prozesses 

J) Beschreibung der Transportprozesse und gegebenenfalls der Abbauprozesse der wirksamen Substanz im 
Organismus, insbesondere bei intravenoser, peroraler oder subkutaner Verabreichung, in dem Rechenpro- 
gramm durch Massentransportgleichungen und kinetische Reaktionsgleichungen, 

K) Verknupfung der nach Schritt C) und D) ausgewahlten Gleichungssysteme und insbesondere numerische 
Berechnung des entstehenden Systems gekoppelter Differentialgleichungen, 

L) Bestimmung des Konzentrations-/Zeitverlaufs der Wirksubstanz in ausgewahlten Kompartimenten. 

2. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, dass der Organismus bzw. die Organismen Sau- 
getiere, aus der Reihe Mensch, Affe, Hund, Schwein, Ratte oder Maus, oder Insekten, insbesondere Raupen, 
oder Pflanzen darstellen. 

3. Verfahren nach Anspruch 1 oder2, dadurch gekennzeichnet, dass der Wirkstoff intravenos, oral, topikal, 
subkutan, intraperitoneal, inhalativ uber die Nase oder Lunge, oder intra-cerebrovaskular (im Fall von Sauge- 
tieren), bzw. intra-haemolymphatisch, topikal, oder oral durch Frafi oder Schlundsonde (im Fall von Insekten) 
appliziert wird, oder die Applikation als Spritzbruhe mit Aufnahme uber Blatter oder Wurzeln (im Fall von Pflan- 
zen) erfolgt. 

4. Ein Verfahren nach den Anspruchen 1 bis 3, dadurch gekennzeichnet, dass die physiologischen Para- 
meter, die die Organismen beschreiben, zeitabhangige Parameter sind. 

5. Verfahren nach den Anspruchen 1 bis 4, dadurch gekennzeichnet, dass die Parameter teilweise einer 
Datenbank entnommen werden, die an das Rechner-System angebunden ist. 

6. Ein Verfahren nach den Anspruchen 1 bis 4, dadurch gekennzeichnet, dass die physiologischen und 
anatomischen Parameter mittels Zufallszahlen im Rahmen einer vorgegebenen statistischen Streuung variiert 
werden (Populationskinetik). 

7. Verfahren nach einem der Anspruche 1 bis 6, dadurch gekennzeichnet, dass fur das Rechenprogramm 
eine hierarchische Model l-Architektur bestehend aus Modulen, die unterschiedliche Funktionalitaten aufwei- 
sen, und 

mindestens zwei Module, die entweder chemische oder biologische Eigenschaften der Substanz beinhalten, 
oder physiologische oder anatomische Informationen uber den Organismus beinhalten, oder Informationen 
uber die Darreichungsform oder die Wirkart beinhalten, verwendet werden. 

8. Verfahren nach einem der Anspruche 1 bis 7, dadurch gekennzeichnet, dass mehrere Modelle von Or- 
ganismen zu einem Gesamtmodell, insbesondere zur Simulation von Drug-Drug-Wechselwirkungen, Mut- 
ter-Foetus-Modell, kombiniertes Insekt-Pflanzen-Modell, zusammengefasst werden. 
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